Experiments on transformation thermodynamics: 
Molding the flow of heat 
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It has recently been shown theoretically that the time-dependent heat 
conduction equation is form- invariant under curvilinear coordinate transfor- 
mations. Thus, in analogy to transformation optics, fictitious transformed 
space can be mapped onto (meta-)materials with spatially inhomogeneous 
and anisotropic heat-conductivity tensors in the laboratory space. On this 
basis, we design, fabricate, and characterize a micro-structured thermal cloak 
that molds the flow of heat around an object in a metal plate. This allows for 
transient protection of the object from heating, while maintaining the same 
downstream heat flow as without object and cloak. 

Transformation optics is a design tool that enables steering the flow of electromagnetic 
waves [U El |3], liquid waves [U |5], elastodynamic waves El El El [TO], or quantum 
mechanical matter waves [Tlj [T2] in an unprecedented manner. In essence, for each of 
these examples, curved space is mapped onto a (meta-)material distribution [T| [T3| 16]. 
For this mapping to be possible, the underlying equations need to be form-invariant 
under curvilinear coordinate transformations. 

Can this basic idea also be translated to other problems such as electrical conduction, 
heat conduction, or particle diffusion? The answer may not appear obvious because 
the time-dependent parabolic differential equations underlying all of these problems are 
distinctly different from the hyperbolic differential equations which govern wave phenom- 
ena. For example, basic phenomena like reflection, scattering, shadows, polarization, or 
interference are crucial aspects in many wave systems, but are completely absent in 
scalar parabolic differential equations. However, conformal maps have already been 
used for more than a century in two-dimensional stationary potential-flow problems in 
hydromechanics to transform boundaries (not materials though) |14j . The case of time- 
independent electric conduction was studied theoretically in 1984 by Kohn and Vogelius 
(following an observation by Tartar) [15] and in 2003 by Greenleaf and coworkers [TO]. 




Figure 1: (a) Square lattice of electrical or thermal resistors, (b) Illustration of the 
curvilinear coordinate transformation underlying the cloak to be shown in 
Fig. [2] The cloak region is shaded in gray with the inner and outer radius R\ 
and i?2 marked. Note that the effective local conductivity varies in both the 
azimuthal and the radial direction. For an observer measuring heat or electrical 
flux outside the shaded region, configurations (a) and (b) are indistinguishable. 
This is the essential idea of thermal cloaking. 

The authors concluded that the two-dimensional tomography problem is not unique. 
This means that one cannot distinguish between, e.g., (i) a homogeneous conductive 
plate and (ii) a plate with a hole and a specially tailored conductivity distribution 
around it by measuring the electric resistance between any numbers of pairs of points 
outside of the central region. Today, this non-uniqueness is commonly referred to as 
cloaking. Cloaking represents a benchmark example for any transformation approach. 
Only quite recently, experiments have been published for steady-state electric cloaking 
|17| . Experiments have also been published [18] within the time-independent limit of 
the corresponding (elliptic) continuity equation for the heat current density (comprising 
a heat shield, but no cloak). However, the more general thermo dynamic problem has 
only been treated theoretically so far [19]. In the dynamic case, the spatial profile of the 
effective specific heat enters explicitly, whereas it is strictly irrelevant in the static case. 
Interestingly, the time-dependent heat-conduction equation is mathematically equivalent 
to the time-dependent particle-diffusion equation. 

In this paper, we experimentally realize a thermodynamic cloak for what we believe 
to be the first time. Movies of the transient temperature profiles in a micro-structured 
copper plate taken with an infrared camera are in good agreement with theory. 

The basic idea of transformation thermodynamics and transformation conduction is 
illustrated in Fig. [I] The resistors connecting adjacent Cartesian grid points in panel (a) 
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can either be Ohmic or thermal. For metals, according to the Wiedemann-Franz law, 
the two types of resistors are simply proportional to each other at a given temperature. 
Upon performing the coordinate transformation given in [I] in panel (b), the grid lines 
become distorted and an empty region is opened in the middle. However, this distortion 
does not influence the flow of the electrical or the heat current outside of the structure, 
i.e., for r > R 2 , at all, provided that adjacent grid points are still connected by identical 
resistors (and assuming, as usual, that the connecting wires have vanishing resistance). 
Thus, the circular hole in the middle can be cloaked. To realize such a fictitious resistor 
distribution by an actual material distribution, it is interesting to note that the area 
density of resistors in panel (b) of Fig. [I] is neither homogeneous nor isotropic. Close to 
the inner edge of the cloak at radius r = R\ in panel (b), the density of parallel resis- 
tors providing current flow in the azimuthal direction is large; hence the electric or the 
thermal azimuthal conductivity is large. In contrast, only few parallel but more serial 
resistors in this region provide conduction in the radial direction. Hence, the radial con- 
ductivity is smaller and even vanishes towards the inner radius r — R±. Mathematically 
|19j . Guenneau and coworkers derived the following approximate form for the azimuthal 
component 



of the heat conductivity tensor k in the radius interval R 2 > r > R\. Here, /^o is the 
isotropic heat conductivity of the surrounding of the cloak. Furthermore, they found 
[T9] that the product of the effective mass density p and the effective specific heat c is 
approximately spatially constant, i.e., pc = const. 

These material distributions still need to be mapped onto a practical metamaterial 
micro-structure exhibiting such pre-described anisotropic heat conduction. Inspired by 
our previous experiments on mechanics [10| , we choose 10 equally thick rings with al- 
ternating large and small effective heat conductivities with a contrast that gradually 
decreases from the inner radius R\ to the outer radius R 2 . The heat conductivities k\ 
of the rings (i = 1-10 from the inside to the outside of the cloak) are determined using 
formulae (24) given in Ref. [19J to obtain the best match to and ^ in the region 
with R 2 >r> Ri. In units of W/(Km), this leads to ^ = 0.15, k 2 = 394.0, ^3 = 2.91, 
k 4 = 390.0, k 5 = 11.26, k 6 = 382.7, k 7 = 19.02, k s = 375.0, k 9 = 26.38, and k 10 = 367.6. 
Intuitively, we exploit the fact that heat conduction is large along the direction of a ring 
with large heat conductivity, i.e., in the azimuthal direction, whereas heat conduction is 
smaller in the perpendicular radial direction. The effective anisotropy depends on the 
contrast within a pair of rings. To act as an effective material, the spacing between 
the rings needs to be small compared to the typical scale of temperature gradients, i.e., 
small compared to the thermal diffusion length. 

To obtain these different heat conductivities, we design a composite structure made 
by drilling holes into a copper plate and filling them with polydimethylsiloxane (PDMS). 




(i) 



and the radial component 
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Figure 2: (a) Blueprint of the designed thermal cloak (only one half is shown due to sym- 
metry). The black regions are bulk copper, the white regions polydimethyl- 
siloxane (PDMS) with heat conductivities of 394W/(Km) and 0.15 W/(Km), 
respectively. As illustrated in Fig. [I](b), the radial and azimuthal components 
of the effective local heat conductivity follow a recipe pre-described by the co- 
ordinate transformation, (b) Photograph of the fabricated thermal cloak and 
(c) scheme of the experiment, results of which are shown in Fig. p) 



4 



The corresponding (room-temperature) bulk heat conductivities are kqu — 394 W/ (Km) 
and ftpDMS — 0.15W/(Km). The latter is more than three orders of magnitude smaller 
than the former. The metal area fraction fa in ring number i is obtained from the 
effective- medium formula K{ — fi^Cu + (1 — fi) ftpDMS- The product pc, which ought 
to be constant, cannot be adjusted independently in our approach. The bulk values 
for copper and PDMS are pc = 3.4MJ/(Km 3 ) and pc = 1.4MJ/(Km 3 ). We find, 
however, that the local effective specific heat merely varies by about 30 % in the radius 
interval [R\,R<2\ for our parameters. To compensate for this approximation, based on 
a numerical optimization, we choose the effective heat conductivity of the surrounding 
= 85W/(Km). 

The absolute size of all mentioned features is not relevant for the behavior of the 
thermodynamic cloak. If we replace r — >> sf and t — > s 2 t with a dimensionless factor 
5, the heat conduction equation remains unchanged. This means that, for example, 
reducing the size of all lateral features by a factor of 10 leads to a thermal time constant 
which is smaller by a factor of 100 — of course provided all material properties remain 
unchanged. We choose R± = 2.5 cm and i?2 = 5 cm, leading to a ring thickness of 2.5 mm. 
The thickness of the plate is not important for the heat-conduction problem of interest, 
but the thickness does influence the strengths of artifacts. For large plate thicknesses, 
the plate's surface-to- volume ratio is large and, hence, heat conduction/convection to air 
is small. However, at some point the overall heat capacity of the plate is so large that 
it becomes difficult to define temperature "baths" , which ought to have heat capacities 
much larger than that of the plate. We find that a plate thickness of 2 mm provides a 
good trade-off. 

In the fabrication of the structure, we start with a bulk copper plate and machine the 
pattern shown in Fig. [2](a) into it. For some of the rings, hexagonal arrays of holes (with 
0.5 mm diameter) are drilled. For other rings, the very low metal filling fractions inhibit 
doing that. In these rings, we rather realize the low metal filling fraction by azimuthally 
evenly distributed thin radial bars. This step makes heat conduction in these individual 
rings anisotropic. However, our experiments show that this approximation has little 
if any detrimental effect on the performance of the system. To compensate for all 
approximations made, we adjust the metal filling fraction in the surrounding. Based 
on an experimental optimization, we choose an hexagonal lattice of holes with lattice 
constant 1.76 mm and diameter 1.5 mm, leading to fo = 34.2%. Thereafter, the holes 
are filled with PDMS, which is polymerized afterwards. Furthermore, we have thermally 
isolated both surfaces of the composite plate by an approximately 100 (xm thin layer of 
PDMS. As a result, heat conduction/convection by air is significantly reduced. 

To measure the transient temperature profiles via Planck thermal emission using a 
conventional infrared heat camera (Flir A320), it is important to achieve nearly 100% 
absorbance = emissivity according to Kirchhoff 's law. Fortuitously, in sharp contrast to 
the highly reflective copper itself, the 100 \im thin PDMS layer on the top of the plate is 
nearly "black" for the wavelengths seen by the thermal heat camera. Thus, we need no 
additional coating of the surface. The absolute temperatures in the temperature profiles 
shown below are derived from the camera data assuming 99 % emissivity (the standard 
value). A photograph of one of the fabricated structures is depicted in Fig. [2](b). Local 
heating on the left-hand side of the plate is achieved by a tank filled with hot water 
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Figure 3: Measured temperature distributions at different times t as indicated. At 
t = s, we start from a homogenous room-temperature profile. Results for 
the complete thermal cloak shown in Fig. [2](b) are shown in the left-hand side 
column, results for a simplified reference structure with only a thermally iso- 
lating ring around the central copper region are depicted in the right-hand side 
column. The white curves are iso-temperature lines (equidistant in steps of 
3° C) corresponding to the temperature profiles depicted on a false-color scale. 
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(see panel (c) in Fig. [2]). The heat capacitance of the filled tank is about two orders 
of magnitude larger than that of the composite plate, qualifying it as a temperature 
bath in the thermodynamic sense. The right-hand side of the plate is immersed in an 
identical tank filled with room temperature water. The remaining two sides of the plate 
are left open (i.e., isolated). This is possible because, unlike waves, the flow of heat is 
not reflected at an open edge. 

Figure [3] shows results of our experiments for the complete cloak (left-hand side col- 
umn) and for a simplified isolating structure (right-hand side column) at various times 
t (different rows) after exposing the plate to the heat baths at t = Os. Prior to that, the 
plate was kept under room temperature conditions for a sufficiently long time such that 
it exhibited a spatially constant temperature at t = Os. The temperature in the plane 
of the plate is shown on a false-color scale. The cloak obviously successfully fulfills its 
task in that the central region is colder than its surrounding. Furthermore, the white 
iso-temperature lines on the downstream side of the cloak in Fig. [3] are nearly verti- 
cal, indicating a temperature distribution as if nothing was there. However, one must 
be cautious because part of the desired and achieved effect can also be obtained by a 
simple thermal isolation around the circular central region. Thus, we have fabricated 
a corresponding reference sample with a solid copper region and the first PDMS ring 
around it (parameters like for the complete cloak). The surrounding is a homogeneous 
perforation just like for the cloak, i.e., this reference structure has no further rings. The 
result depicted in the right-hand side column of Fig. [3] also shows a central region that 
is cooler than its surrounding. This finding is very much different from light waves or 
mechanical waves due to the fundamental difference between hyperbolic and parabolic 
differential equations emphasized in the introduction. The optical counterpart of the 
weakly heat conducting PDMS ring would be an opaque metal film. In mechanics, the 
analog would be a rigid wall. As a result, one gets scattering, shadowing, and wave- front 
distortions. These effects have, e.g., been observed in Ref. [10] for a comparable two- 
dimensional scenery. In contrast, in Fig. |3| all of these effects are absent. One merely 
finds distortions of the iso-temperature curves, which may be seen as the counterpart of 
wave fronts. The iso-temperature curves in the homogeneously perforated region near 
the center of the reference structure are curved towards the center. Here, the heat flux 
is simply perpendicular to the iso-curves. In sharp contrast, the iso-temperature curves 
for the complete cloak in Fig. [3] exhibit the opposite curvature, namely away from the 
cloak, both on its left (leading to an apparent negative diffusivity at each time step) 
and on its immediate right. Here, the heat flux cannot easily be seen from the images 
because it results from the product of the local anisotropic heat conductivity tensor and 
the local temperature gradient. 

Generally, our experiments in Fig. [3] agree well with the corresponding numerical 
calculations (Comsol) shown in Fig. Hi The heat baths are mimicked by fixed tem- 
peratures. Furthermore, the calculations assume isotropic homogeneous rings with the 
prescribed heat conductivities given above and do not account for heat conduction to 
the surrounding air. The latter aspect likely explains the remaining discrepancies. We 
refrain from showing the modeling accounting for this aspect because it is not part of 
the cloak design. 

In conclusion, using the theoretical concepts of transformation thermodynamics, we 
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Figure 4: Calculated temperature distributions shown as the experimental results in 
Fig. [3j For the cloak in the left-hand side column, the ring structure is in- 
dicated for clarity. 
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have realized a cloak that molds the dynamic flow of heat around an object as if no object 
was there. Such thermal cloaks might find, e.g., applications in temporarily protecting 
sensitive regions in an electrical circuit or chip from excessive heating. More importantly, 
our experiments also demonstrate that the ideas of transformation optics go well beyond 
systems exhibiting waves and also work for heat conduction. Realizing inhomogeneous 
and strongly anisotropic thermal metamaterials is far easier than in other systems. In 
the present paper, the material contrast (i.e., in heat conductivity) of the constituent 
materials is larger than 2600. For comparison, in our previous mechanical experiments 
[TO] , the material contrast (i.e., in phase velocity) was about 39. For visible- frequency 
off-resonant dielectrics in optics, this ratio is below 3. 

We thank Frank Landhaufier, Mario Nusche, and Johann Westhauser (KIT) for help 
in the fabrication, Christoph Kaiser (KIT) for lending to us the infrared heat camera, 
and Jonathan Miiller (KIT) for stimulating discussions regarding the cloak design. We 
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